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Semi-leptonic and leptonic decays of B-mesons are important probes for testing SM and theories 
beyond it because of their relative cleanliness and far less theoretical uncertainties. In semi-leptonic 
decays based on quark level transition b — > s£ + £~ apart from branching ratio one can study many 
other (possible) observables associated with final state leptons like, lepton pair Forward Backward 
asymmetry, lepton polarization asymmetries etc. But as proposed recently if we can tag the B- 
meson than one can measure the polarization asymmetries of both the leptons. Here we will study 
the polarization asymmetries of both the final state leptons in SM and Minimal Supersymmetric 
Extension (MSSM) to it 



PACS numbers: 



I. INTRODUCTION 



The rare B-meson decays induced by Flavor Chang- 
ing Neutral Current (FCNC) b — > s(d) transitions arises 
only at loop levels in Standard Model (SM). Because 
these transitions occur at loop level, hence provides useful 
tests of the detailed structure of the theory at the levels 
where GIM (Glashow-Iliopolus-Maini) cancellations be- 
comes very important. Also in most of the extensions 
of SM, loop graphs with new particles (most of the ex- 
tensions of SM predict existence of some new particles 
like SUSY predicts whole set of SUSY particles) can con- 
tribute to the same order (in case of pure dileptonic decay 
of B s meson SUSY can give a contribution which is many 
orders greater than SM contribution [l|,|2(). ^ n particular 
the process B — > X s ^ and B — * X s £ + £~ are experimen- 
tally very clean and can possibly be more sensitive to 
any new Physics beyond the SM. The new physics effects 
in the rare decays can come in two ways, one via the 
modifications of the existing (in SM) Wilson coefficients 
and other via the introduction of some new operators 
(accompanied by new coefficients) Q- Ql • The inclusive 
decay mode (B — > X s £ + £~ ) can be more sensitive than 
the radiative decay mode (B — > X s j) for testing any 
new physics model because in inclusive decay mode many 
more kinematical distributions like, lepton pair forward 
backward (FB) asymmetry, leptons polarization asym- 
metry etc. can be measured. 

Various kinematical distributions of the inclusive mode 
have been studied in many earlier works 0, . As also 
been proved in many of the works that lepton polariza- 
tion asymmetry of final state leptons can give us useful 
information to fit parameters of SM and constraint new 
physics models 0, 0. But recently it has been noted 
down, in case of inclusive decays by Bensalem et.al, 
that one can in principle observe many more observables, 
like the double polarization asymmetries (polarization 
asymmetries when both the leptons are polarized), which 
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would be useful in further testing of the SM and probing 
physics beyond it. They inferred that even if we won't be 
able to tag b but can observe the final polarization state 
of both the leptons than also we can have more distribu- 
tions than what we have if we would only be measuring 
the single lepton polarization asymmetries. But if we can 
also have b-tagging along with the measurement of the 
polarization of both the final state leptons, then many 
more kinematical distributions would be available to us, 
which could be useful in testing the structure of effective 
Hamiltonian and hence the physics underlying it. In our 
earlier work we tried to estimate as to what would 
be the form of various double polarization distributions 
within SM for the exclusive mode B — > K*t + t~ and 
how these distributions gets modified by switching on the 
SUSY effects. This work is the extension of the work of 
Bensalem et.al., 0. In their work they restricted them- 
selves to SM operators. In this work we will extend the 
SM operator basis and will try to observe the aftereffects 
of this on the various double polarization asymmetries. 

The quark level transition we would be interested in 
the present work is b — > s£ + £~ . But as has been very 
well emphasized in literature that if we consider the Su- 
persymmetric (SUSY) extension of the SM, than we have 
to extend the SM list of operators to include the oper- 
ators arising from the exchange Neutral Higgs Bosons 
(NHBs) 0,0,0, IE I 11 case °f pure-dileptonic decays 
of B-meson (B s -> £ + £~) the NHBs can change the SM 
predictions by many orders of magnitude U El- As 
the couplings of the NHBs to leptons is proportional to 
the lepton mass, hence these effects becomes more impor- 
tant if t = fj,, t. The SUSY extension of the SM predicts 
existence of two new operators, which were not present 
within SM. These new operators are responsible for the 
orders enhancement of the branching ratio of pure dilep- 
tonic decay mode of B-meson. In our analysis we will 
be going to consider these operators and will try to esti- 
mate the dependence of the various double polarization 
asymmetries on MSSM parameters. 

The paper is organized as follows : In section [H] we 
will present the effective Hamiltonian we are consider- 
ing. We will then write the matrix element for the quark 
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FIG. 1: Branching ratio of B — > X,t t with scaled in- 
variant mass of dileptons. Parameters of mSUGRA model 
are : m = 400GeV, M = 500GeV, ,4 = 0, tan/3 = 40 and 
sgn(n) is positive. The additional parameter for rSUGRA is 
m A = 270GeV 



level process b — > s£ + £~ . In section ITTT1 we will give the 
definition and the analytical results of the double polar- 
ization asymmetries. Finally, we will conclude in section 
II VI bv discussing our numerical results. 



II. 



EFFECTIVE HAMILTONIAN 



We follow the convention followed in 0, || to write 
down the matrix element and invariant mass spectrum. 
The inclusive decay mode (B — > X s £ + £~ ) is modeled by 
partonic process b(pb) — > s(p s ) + £ + {p + ) + £~{pJ). This 
sort of modeling is actually the leading order calculation 
in X/mb expansion. By integrating out the heavy degrees 
of freedom from the full theory (MSSM here) we can 
get the effective Hamiltonian describing the semi-leptonic 
decay s£+£~ [IS 0,0 : 
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V tb V* \ >:C i O i + ^C Qi Q i 

i=l 



(2.1) 



where Oi are current-current (i=l,2), penguin (i = 
1,...,6), magnetic penguin (i=7,8) and semi-leptonic (i 
= 9,10) operators and Ci are the corresponding Wilsons. 
They have been given in 0, ^| . The additional opera- 
tors Qi (i = 1,. . . ,10) and their Wilson coefficients JCn i ) 
which arises due to MSSM diagrams are given in 

Neglecting the mass of s-quark, the effective Hamilto- 
nian gives the matrix element : 
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olGf 
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+Cl ff (s^P L b)(£^£) + C w (s 7fl P L b)(i^ 75 i) 
+C Ql (sP R b)(££) + C Q2 (sP R b)(£ l5 £) 1 (2.2) 
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FIG. 2: Total decay rate with tan/3 in mSUGRA model. 
Other parameters m = 400GeV, A = 



where q is the momentum transfer to the lepton pair 
given as q = p + +p_ , where p_ and p + are the momentas 
of t~ and £ + respectively. VtbV t * s are the CKM factors 
and P LM = (1 T 7s)/2 • 

The matrix element A4(b — » s£~£ + ) is the free quark 
decay amplitude but it has some long distance ef- 
fects also which comes due to the four-quark operators 
(£ + £~s\Oi\b) where i — 1,...,6. These effects are ab- 
sorbed into the redefinition of short distance Wilson co- 
efficients according to the prescription given in many ear- 
lier works 0, El • The long distance effects are because 
of the the cc resonance contributions. These resonances 
can be taken into account by using Breit-Wigner ansatz 
H H3 by which we add a term to Cg • The resonance 
contribution to Cg is : 
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where the symbols are explained in Kruger & Sehgal |(| • 
The phenomenological factor k has to be introduced to 
reproduce the correct branching ratio B{B — > J/^>X S ) — > 
X s £ + £~). We will be going to take this factor to be 2.3 
for our numerical analysis. 

From the expression of the matrix element given in 
ean. (|2.2l) we can get the dilepton invariant mass distri- 
bution as : 
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FIG. 3: Total decay rate with m A GeV in rSUGRA model. 
Other parameters m = 400GeV, M = 500 GeV, A = 
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FIG. 5: Vln with s with all the parameters same as in Fig0 




FIG. 4: Vll with s with all the parameters same as in Fig. 
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Now ready with the expression of the invariant mass 
spectrum (including the scalar exchange effects) we will 
analyze various double polarization asymmetries in next 
section. As has been emphasized in many texts that 
the r-polarization asymmetries in inclusive decay B — > 
X s t + t~ can be a very useful probe of the structure of 
effective Hamiltonian and hence the underlying theory. 
So for our analysis we will consider the inclusive decay 
channel with r leptons in final state. 



III. 



LEPTON POLARIZATION ASYMMETRIES 



the convention as given in many earlier works 00113 ■ To 
evaluate the polarized decay rates we have to introduce a 
spin projection operator defined by N = 1/2(1 + j5/6 x ) 
for t~ and M = 1/2(1 +75 for £+ where x = L,N,T 
and corresponds to the longitudinal, normal and trans- 
verse polarization asymmetries respectively. Firstly we 
define the orthogonal unit vectors S x for £~ and W x for 
t + in rest frames of £~ and £ + respectively as : 
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where P-,p+ and p s are the three momentas of £~,£ + 
and strange (s) quark in center of mass frame (CM) frame 
of £~~£ + respectively. Now from the rest frame of respec- 
tive leptons we boost the four vectors S x and W x to the 
dilepton CM frame. Only the longitudinal vectors, Sl 
and Wl, will get boosted by the Lorentz transformation 
to CM frame of £~£ + to a value : 
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In this section we will evaluate the double lepton polar- 
ization asymmetries, i.e. where polarization of both the 
leptons is begin measured. For this we have to define the 
polarization vectors of l~ and i + . We will be going to use 



where Eg is the energy of any of the leptons (both have 
same energy) in dileptonic CM frame. 

Now we can define the double polarization asymme- 
tries as : 
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where the sub-index x, y are L,N or T. 

We can get the expressions of double polarization 



asymmetries as : 
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where A is given in ean. H2.5|l . 



IV. NUMERICAL ANALYSIS, RESULTS AND 
DISCUSSION 



In this section we will be going to discuss our numerical 
analysis and the results of our numerical analysis. Firstly 
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FIG. 6: Vlt with s with all the parameters same as in FigQ 
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FIG. 7: Vnn with s with all the parameters same as in FigQ 



we will be going to list the SM values of the branching 
ratio of B — ► X s t + t~ and various double polarization 
asymmetries we have presented in previous sections in 
Tabled 

We have analyzed the SUSY effects on the various ob- 
servables which we have listed in previous section. MSSM 
is although the simplest and the one having least number 
of parameters, extension of the SM. But still it has large 
number of parameters which makes it difficult to do phe- 
nomenology with it. But we do have some models, like 
Dilaton, moduli, mSUGRA, . . . etc, which reduces this 
vast parameter space to a manageable level. In our nu- 
merical work we will be going to use one of the more pop- 
ular unified model called Supergravity (SUGRA) model. 
The main feature of SUGRA (we will call this as minimal 
SUGRA or mSUGRA) model is that a unification of all 
the scalars, fermions and coupling constants is assumed 
at GUT scale. So effectively in mSUGRA framework we 
have five parameters which are : m (unified mass of all 
the scalars) , M (unified mass of all the gauginos) , A (uni- 
versal trilinear coupling constant), tan/3 (ratio of vev of 




FIG. 8: Vtl with s with all the parameters same as in Fig0 




FIG. 9: Vtt with s with all the parameters same as in FigQ 



the two Higgs doublets) and sgn^n) 1 . 

But as we know that universality of scalar masses at 
GUT scale is not a necessary condition of SUGRA [T^| . 
One can in principle have a different unified mass of the 
squark sector and Higgs sector at GUT scale. We will 
explore this scenario also, which we will call as relaxed 
SUGRA (rSUGRA) model. The additional parameter 
which we will be going to have in this model we will take 
as mass of pseudo-scalar Higgs boson (m^)- 

We will be going to work in the high tan/3 region of the 
SUGRA parameter space because its only in this region 
the effects of NHBs becomes more prominent. There are 
some constraints on the MSSM parameter space from ex- 
perimental observations of B — ► X s ~/ 13]. In our analysis 
we will be going to consider only that region of parameter 



1 the convention about sgn(fi) we will be going to follow is that it 
appears in chargino mass matrix with a +ve sign 
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Br(B -> X s r + r" ) 




Vln 


Vlt 


Vnl 


Vnn 


Vnt 


Vtl 


Vtn 


Vtt 


1.29 x 10~ 7 


- 0.082 


- 0.137 


- 0.136 


0.136 


0.122 


- 0.017 


- 0.397 


0.017 


-0.216 



TABLE I: Standard Model predictions of the observables for B — ► X s r + r 



space which satisfy the 95% CL bound : 

2 x 1CT 4 < Br(B -► Y s7 ) < 4.5 x 10~ 4 

which is agreement with CLEO and ALEPH results. 

We will also present the numerical results of the aver- 
age polarization asymmetries. The averaging procedure 
which we will be going to use is defined as : 



(V) 



(m B -m K *) 2 /m 2 B 

T—ds 

(3.646+0.02) 2 /m| « s 



(m B -m K * ) jm B 
(3.646+0.02) 2 /m^ ds 



(4.1) 



-ds 



which means in integrating the observables we will be 
going to consider the region of dilepton invariant mass 
which is above the first charmonium resonance (after the 
threshold of the process b — ► st + t~). We will be not be 
going to present the result of Vnt (and Vtn) as they are 
very small. 

We have presented the plots, of the various kinemati- 
cal variables, we presented in sections [H] and IIIII with 
the (scaled) invariant mass of the dilepton. In FigQ] 
we have shown the variation of the branching ratio of 
B — > X s t + t~ with dilepton invariant mass in three dif- 
ferent models we have considered, namely SM, mSUGRA 
and rSUGRA. In Fig. |2 il and we have 
plotted the distributions of various double polarization 
asymmetries with dileptonic invariant mass. As we can 
see from these figures that SUSY can change these dis- 
tributions substantially over the whole kinematically al- 
lowed region. We have done the detailed scanning of 
the mSUGRA and rSUGRA parameter space and have 
presented the results of the integrated polarization asym- 
metries, as defined in eqn.gT}. In Fig. ES El OH! HHI 
and|2D]we have shown the plots of various integrated po- 
larization asymmetries with tanfi for various values of M 
(the unified gaugino mass at GUT scale) in the mSUGRA 
model. As we can see that there is a substantial changes 
in the SM and mSUGRA model prediction. Specially for 
Vll where depending upon the mSUGRA parameters it 
can even change its sign. In Fig. [l^fT^lfTBlflTlfT^l aridl^ll 
we have plotted various integrated polarization asymme- 
tries as function ttia (mass of pseudo-scalar Higgs boson) 
for various values of tan (3. There also one can observe 
major changes in the predictions of various integrated 
polarization asymmetries as compared to SM values. 
As has already been mentioned in many works 0, 0, 
] that polarization asymmetries are useful in find- 
ing out the structure of the effective Hamiltonian and 
hence the physics underlying it. As also been argued 



by Bensalam et al. 0, if we work within the SM then 
for inclusive decay B — * X s £ + £~ we have five theoretical 
parameters which are : four Wilsons (CV, C10, real and 
imaginary part of Cg ) and mj,, they can in principle be 
completely determined using the three r~ polarization 
asymmetries, the total decay rate and the FB asymmetry 
i.e. five observables. But as known from earlier results 
0, |(| the normal polarization asymmetry is very small 
and hence we have to consider the polarization asymme- 
try of r + also. Also the measurement of FB asymmetry 
requires b-tagging. So if we have a untagged sample then 
the observables which we have are : decay rate, two polar- 
ization asymmetries of t~ (we are neglecting the normal 
one being very small) and two polarization asymmetries 
of t + (again neglecting the normal one). So we have five 
parameters and five observables which should in princi- 
ple be sufficient. But there is no other constraint which 
will give us the cross-check of SM. So along these argu- 
ments Bensalam et.al. constructed double polarization 
asymmetries which give us large number of observables 
to have a good cross-checking of SM. 

Let's now examine the situation if we believe in SUSY 
extension the SM. As has been argued in many works 
01 IMS that if we consider the decay channel b — > st + t~ 
in SUSY extension of SM 2 , then we have to extend the 
SM list of operators by introduction of two new opera- 
tors Cq 1 and Cq 2 . So now our theoretical parameters 
would be seven (five SM one and two new Wilson coeffi- 
cients). So we require more number of observables to fix 
up this new structure of effective Hamiltonian and hence 
the double polarization asymmetries could possibly be 
very useful. 

The model which we are considering is MSSM. In 
MSSM if we assume the MSSM parameters to be real 
then the only source of CP violation is the CKM matrix. 
The transition b — > s£ + £~ is a CP conserving process. 
Let's try to analyze the number of additional observables 
which we now have 3 . For this we consider two possi- 
bilities : one where we can have b-tagging and the other 
where there is no b-tagging. If we say that the double 
polarization asymmetries for B — > X s t + t~ are Vij then 
for the conjugated process B — > X s t + t~ we denote them 
with V^. But in general Vij = ±Vij 4 . 



2 in fact this is true even if we consider the two Higgs Doublet 
model (2HDM) extension of SM 

3 by additional observables we mean the observables due to mea- 
surement of double polarization asymmetries 

4 in fact except for Vnl and Vnt the sign is always positive. For 
Vnl an d Vnt the sign is negative 
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So if we consider the first case where there is no b- 
tagging. In a untagged sample we can measure four ad- 
ditional observables namely Vll, VnNi Vtt and (Plt + 
Vtl)- 

But if we consider that we can tag b in that case we 
have all the nine double polarization asymmetries avail- 
able to us which will give us very useful probes into the 
structure of effective Hamiltonian. 
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APPENDIX A: INPUT PARAMETERS 

m B = 5.26 GeV , m b = 4.8 GeV , m c = 1.4 GeV 
m r = 1.77 GeV , m w = 80.4 GeV , 
m z = 91.19 GeV , V tb V t * = 0.0385 , 
a = m ' T B = 4.22 x 10- 13 GeV 
G F = 1.17 x 10" 5 GeV" 2 
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FIG. 10: Integrated Vll with tan(5 in mSUGRA model. Parameters same as given in Fig|2l 
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FIG. 11: Integrated Vll with m^GeY in rSUGRA model. Parameters same as given in Fig. 



SM 

M = 450 
M = 550 
M = 650 



I 

30 



35 40 45 50 



FIG. 12: Integrated Vln with tanf3 in mSUGRA model. Parameters same as given in FigEl 
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FIG. 13: Integrated Vln with m^GeY in rSUGRA model. Parameters same as given in Figl3"l 
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FIG. 14: Integrated Vlt with tan/3 in mSUGRA model. Parameters same as given in Fig. [5] 
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FIG. 15: Integrated Vlt with m^GeY in rSUGRA model. Parameters same as given in Fig. [3 
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FIG. 16: Integrated Vnn with tan/3 in mSUGRA model. Parameters same as given in Fig. [5] 
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FIG. 17: Integrated Vnn with niAGeV in rSUGRA model. Parameters same as given in Fig. |3 
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FIG. 18: Integrated Vtl with tan/3 in mSUGRA model. Parameters same as given in Fig. [5] 
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FIG. 19: Integrated Vtl with mAGeV in rSUGRA model. Parameters same as given in Fig. 
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FIG. 20: Integrated Vtt with tanfi in mSUGRA model. Parameters same as given in Fig. [2] 
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FIG. 21: Integrated Vtt with m^GeV in rSUGRA model. Parameters same as given in Fig. 



